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Abstract 

Two-proton correlations at small relative momentum q were studied in the eA( 3 He, 4 He, C, Fe)^ 
e'ppX reaction at E$ = 4.46 GeV using the CLAS detector at Jefferson Lab. The enhancement 
of the correlation function at small q was found to be in accordance with theoretical expectation. 
Emission region sizes were extracted and proved to be dependent on A and proton momentum. 
The size of the two-proton emission region on the lightest possible nucleus, He, was measured for 
the first time. 

PACS numbers: 21.65.+f, 25.10,+s 
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One of the outstanding issues in nuclear physics is the nature of dense and(or) hot nuclear 
matter . There are strong experimental indications |3j, |4| that density fluctuations of 
nuclear matter manifest themselves in so called "cumulative processes" in which particles 
are produced in the kinematical region forbidden for interactions with a single motionless 
nucleon, and hence more than one target nucleon must be involved. Cumulative particle 
spectra remain unexplained when finite temperature Fermi-gas momentum distributions 
are considered leading to the association of the reaction strength in this kinematic 
region with density fluctuations or correlations. These objects can be described in various 
wa y s but a „ autho, coas.de, t hem to be flue—. lu this paper, we wfl, 

not rely on a specific model, and, following Blokhintsev Efremov jg] and others, will, for 
simplicity, refer to this type of object as a "flucton". The production of an energetic nucleon 
pair from a nucleus is also an example of a cumulative process, and therefore can be used 
to study fluctons. 

Pairs of nucleons can also be produced due to rescattering of the emitted particles on 
other nucleons in the nucleus. Cascade calculations also fail to describe the whole set of 
experimental data, but rescattering can affect experimental spectra and particle correlations. 
The relative importance of rescattering processes depends on the mass number A of the 
nucleus. We believe that extrapolation to the smallest A will provide reliable information 
on the true properties of the flucton. 

To estimate the density of the flucton, one needs to measure its size and the number of 
contributing nucleons, the minimum number of which can be determined from the kinemat- 
,e, The fluctou s.e is expected to be co_ate wu, the s iz e of a uuOeou BflS. 

The two-particle correlations at small relative momentum q = p[ — p 2 {p[ and p~2 are 
;he individual proton momenta in the pair reference frame) are sensitive to the source size 
14J, [15[ (see also the reviews |16|). We will use the term "femtoscopy" (1 fm = 10 15 
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m) for the study of source sizes within nuclei in analogy with microscopy. 

The two-proton correlations at small q was theoretically described in jl^l ll^ . The inter- 
ference of identical particles jl^ . 3] , as well as Coulomb and strong final state interactions 
(FSIs) were taken into account. Strong FSIs are dominant, causing the increasing of 
the pair production cross section near q ~ 0.04 GeV/c. The intensity of the effect depends 
inversely on the root mean square radius trms of the source from which the protons are 
emitted. 
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It should be noted that here we understand the FSI as the interaction in the two-proton 
system at small relative momentum only. The interaction time in this system is much larger 
than the characteristic collision time and so this system can be considered in isolation of 
other particles and described by the same wave function as in the scattering problem (up to 
the opposite direction of the relative momentum vector). For proton interactions with other 
particles during the collision process we use another term - rescattering. Rescatterings are 
generally characterized by much higher momentum transfers and correspondingly shorter 
time scales. They are essentially localized and can be considered as new emission points. 
FSIs are our "tool" to measure the ffucton size while the rescatterings wash out the original 
emission region and thus distort this measurement. 

Although femtoscopy has been used widely to study a number of processes 
{hh,e + e~,AA (rJ) this is not the case for the cumulative process. Hadroproduction data 
exist for carbon and heavier nuclei but lepton-nucleus data are scarce in any kinemat- 

ical domain 0,0]. In 0] the size of the pion emission region was studied in high energy 
z/D interactions. In [20] data on two-proton and two-pion correlations were obtained in e 16 
interactions at 5 GeV. The scattered electron was not identified; the data correspond to a 
small Q 2 value and to v ~ 1.5 GeV. The measured source size proved to be commensurate 
with the nuclear size and showed a tendency to decrease with particle momenta. 

We present here our study of the correlation between two detected protons with small 
relative momenta in eA( 3 He, 4 He, 12 C, 56 Fe)-^ e'ppX reactions, for an incident electron 
energy of 4.46 GeV. The measurements were performed with the CEBAF Large Acceptance 
Spectrometer (CLAS) 21| in Hall B at the Thomas Jefferson National Accelerator Facility. 
The CLAS detector is a six-sector toroidal magnetic spectrometer. The magnetic field is gen- 
erated by iron-free superconducting coils. The detection systems consist of drift chambers 
to determine the trajectories of charged particles 0], scintillator counters for time-of-flight 
measurements 0], Cherenkov counters to distinguish between electrons and negative pi- 



ons 



24j . and electromagnetic shower calorimeter to identify electrons and neutrons [21 



The CLAS was triggered on scattered electrons detected in the calorimeter with an energy 
above 1 GeV. 

Run conditions are described in detail in Ref. 0]. Events with transferred energy v be- 
tween 0.5 and 3.5 GeV and transferred 4-momentum squared Q 2 between 0.6 and 5 (GeV/c) 2 
were selected for analysis. Protons in the momentum range from 0.3 to 1.0 GeV/c were se- 
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lected for the analysis. The angle 9 between the direction of the virtual photon 7 and the 
direction of the detected proton was from 0° to 115°. The analysis was performed for events 



with at least two detected protons. Misidentification o 
In this article we shall use the "mixing" procedure 
calculations, i.e. 



electrons or protons was negligible. 
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for the correlation function(CF) 



N r (q,p) 

R[q ' p) = Nj^fy (1) 

where q — \q\, p = \p\ and p = {p\ + j>2)/2; N r and N m are the numbers of proton pairs from 
the real events and those combined from protons taken from different events, respectively. 
Secondary particles are boosted in the direction of the virtual photon momentum. We select 
the mixed-pair protons from events for which the magnitude of the momentum difference 
of the scattered electrons \p~li — p^\ is less than q . We studied the dependence of the N m 
distribution on the value of go, and found this dependence negligible for q <0.2 GeV/c; a 
cut q Q <0.2 GeV/c was applied for the N m distributions. Pairs of tracks hitting a single 
scintillator were not included in our analysis, because they have ambiguous time-of-flight 
values. 

The ability to detect two tracks with small relative momentum is limited because both 
particles hit the same or neighboring detector cells. As a rule, the probability for losing at 
least one of the two tracks is higher for close tracks. A detailed study of the close-track 
efficiency e(q) in CLAS has been done in Ref. It depends on track curvature and then, 
for a fixed nominal magnetic field, on the proton momentum and emission angle in the 
laboratory system. The dependence of e(q) for the mean momentum and emission angle is 
shown in the insert of Fig. ^ 

Fig. n shows R{q) for the 3 He, 4 He, and Fe data corrected for close-track efficiency s(q), 
"long-range" correlations(LRC), and momentum resolution. The data are averaged over 
proton momenta. LRCs arise mainly from momentum conservation for the real events which 
is not a requirement for mixed pairs. It results in a smooth increase of R with q, which 
reflects the fact that due to momentum conservation the probability of particles emission in 
the same direction is smaller than that in the opposite direction. Empirically, LRC can be 



parameterized 
is a parameter 



yyR oc exp(6cos , 0), in which ip is the angle between the two protons and a b 
281 ] . The parameter b for different A and proton pair momenta, was obtained 
from a fit to the data in which the region of the effect at small q < 0.2 GeV/c was cut out. 
The correction for LRCs were made by introducing a weight w = exp(bcosip) for mixed 



pairs to reproduce LRC in the N m distribution. The proton momentum resolution within 
the selected kinematic range is estimated to be 5p/p ~ 2 %. Since 6p typically smaller 
than the width of the effects under study measured correlation functions are only slightly 
smeared out by momentum resolution. The momentum resolution corrections were made by 
applying the smearing procedure n times to the measured CF and then by the extrapolation 
of the results to n = —1. 

Fig. [T] also shows the theoretical dependencies of R(q) |3 f° r r rms = 1.6 and 3.0 fm 
calculated within the model of independent one-particle sources taking into account quan- 
tum statistics and FSIs in the two-proton system. The theoretical correlation function is 
then calculated as a square of the wave function (corresponding to the scattering problem) 
averaged over the relative distances of the emitters in the pair rest frame. We assume a 
Gaussian distribution of the emission coordinates in the nucleus rest frame characterized by 
a dispersion = r^ ms /3. We neglect here the emission duration which enters in the longitu- 
dinal component of the relative distance vector through the Lorentz transformation to the 
pair rest frame (the duration is thus effectively absorbed by the parameter r rms ). Both the 
correlation functions and the theoretical curves are normalized to unity for 0.17< q <0.35 
GeV/c. The theoretical approach predicts that the enhancement of R at small q 

is inversely related to the measured size parameter. The peak at q ~ 0.04 GeV/c results 
mainly from the interplay between the attractive s-wave strong final-state interaction and 
the Coulomb repulsion. We compared the results of the calculation of the theoretical CF 

14. llfil l^: i.e. spherical wave approximation (scat- 



for different proton-proton potentials 



tered wave ~ V r )> simple square well potential, Reid [30], and Tabakin [31]. For large 
r RMS values, the correlation function is mainly determined by the solution of the scattering 
problem outside the range of the strong interaction potential, and is therefore independent 
of the actual form of the potential, provided that it correctly reproduces the scattering am- 
phtndes 0Q. 0ur resnhs start to depend on the potent,, choice for ^rms <2 fm. At 
'"rms <2 fm the calculated curves for different potentials also look similar, but the best value 
of trms depends on the version of the potential. In the present work final results for trms 
are presented for the realistic Reid potential, and the difference between results calculated 
for the potential with core (Reid) |3(| and without core (Tabakin) j3l| (~ 3% in r RMS for 
the He data) is taken as the theoretical uncertainty. 

The curves in Fig. ^ represent the best fit to the data by the theoretical curves (the 
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FIG. 1: The two-proton correlation function R for 3 He, 4 He and Fe nuclei. Curves are calculated 
for vrms =1-6 fm (He) and trms =3.0 fm (Fe). The insertion shows the close-track efficiency e(q) 
and its uncertainty beje. 

difference between 3 He and 4 He is negligible) with vrms as a free parameter. The fits in 
Fig. are quite reasonable. This is an indication that the theoretical approach is applicable 
down to a measured size of the order 1.5 fm. The dependencies of R on q for 3 He and 4 He 
(and the best value for t^ms) are the same within errors; the enhancement of R at small q 
for Fe is much smaller. This means trms is larger for Fe than for He. The results for carbon 
(not shown in the figure) lie between He and Fe. 

Experimental systematic errors on r rms arise from the close-track efficiency correction 
(~ 2%), the correction for "long-range" correlations (~ 2%), and the correction for momen- 
tum resolution 1%). Non-identified A particles that decay into pir provide a potential 
background for the measured CF. The cross section for A production is estimated to be 
smaller than 1% of the proton production cross section in the corresponding kinematical 
region. (Additional mass S m ~ rrix + rriA — m p must be produced, and in a cumulative pro- 
cess the cross sections falls exponentially with 5 m , having a slope parameter on the order of 
the pion mass). Therefore, the background from non-identified A — > pir decay is negligible. 
The total systematic experimental errors on trms is about 3%. In the figure statistical and 
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FIG. 2: The size parameter Trms as a function of the mean pair momentum p = \p\ + p2\/2. 
Data [2(J] , which correspond to e 16 interactions at initial energy 5 GeV and Q 2 < 0.1(GeV/c) 2 , 
are shown for comparison. 

systematic errors have been added in quadrature. 

The dependence of trms on p = \p\ + p*2 1 / 2 for different nuclei is shown in Fig. 121 The 
data are averaged over emission angles, statistical and systematic errors have been added 
in quadrature. For 3 He the momentum dependence looks flat, while for carbon and iron it 
decreases with increasing pair momentum. Our results for carbon are in good correspondence 
with the data j^J for electron-oxygen interactions. The values of trms approach the size of 
the nucleus for the lowest value of pair momentum, which seems to be due to the rescattering 
of protons in nuclear matter. The importance of rescattering decreases with proton momenta 
in the chosen momentum range due at least in part to the decrease in the NN cross section. 

We estimate the size of the flucton rj under the assumption that both the primordial 
source size r f and its modification due to rescattering processes contribute to the measured 
size. In the case of helium, the probability of rescattering is much smaller than in heavy 
nuclei. The extracted trms values in 3 He and 4 He are about the same, which is additional 
evidence that rescattering does not affect the helium data within the errors (~ 0.1 fm). 
Therefore, r^Ms in helium 1.6 fm) is an upper estimate of r/. 
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To take into account the possible influence of the rescattering process for helium, we can 
extrapolate the measured sizes as a function of A to the minimum possible target mass, where 
rescattering is not possible. This will provide a lower estimate of Tf, because rescattering 
can only increase the measured size. The minimal target mass (in nucleon mass units) for 
the electro- production of protons(the so-called cumulative number X$ 32]) is determined 
by the kinematics of the process e + X$ ■ m p — > e' + p + m c and is given by: 

_ g + E p -P p cose p7V /l + Q7^ 
(1 - T p /v)m p 

in which E p , P p , m p and T p are the full energy, momentum, mass and kinetic energy of 
the proton, and 6 P1 is the angle between the proton momentum and the virtual photon 
momentum and m c is determined by conservation laws for quantum numbers, baryon number 
in our case. In the limit of large v, X$ approaches the sum of the Bjorken variable Xsj = 
Q 2 /2m p v and the light cone variable a = (E p — P p cos6 P1 )/m p . For a di-proton (a proton 
pair at small relative momentum) the electro-production cumulative number is given by 
Eq. |21 in which E p , P p , T p and 9 P1 now refer to the pair. 

Cumulative production is defined to occur when Xs is larger then unity. Half of our 
proton pairs are produced with Xs > 2; the remaining events are still close to the kinematic 
boundary in the reaction when the mass of the target is the two-nucleon mass. An extrapo- 
lation of the measured sizes to A ~ Xs provides the pair momentum average value 1.2 ±0.1 
fm, where the error arise mainly from the dependence of the result on the extrapolation 
law. The measured CF and then the extracted r RM s could be affected by background from 
the decay of short-lived resonances like the A. Since the proton velocity in the A decay 
reference frame is small v ~ 0.2 and the lifetime is of the order cr ~ 2 fm, this background 



could contribute \J (trms) 2 + ( VT ) 2 ) — r RMS to the measured size, which is less than 0.1 fm. 
Given the maximum possible value of this background the lower estimate for the flucton size 
is 1 fm. Therefore, we estimate the flucton size as rf = 1.3 ± 0.3 fm, which is an average 
of the 1 fm lower estimate and the measured value for He of 1.6 above. It should also be 
noted that the flucton size estimate in 0, Q| was indirect, rather imprecise, and based on the 
model for fitting inclusive data only. This work presents direct measurement of the flucton 
size. 

In summary, the correlations between protons produced in eA interactions at 4.46 GeV 
have been investigated. The data clearly show a narrow structure in the correlation function 
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in the region of small relative momenta (q <0.1 GeV/c) with a peak at q ~0.04 GeV/c which 
is in accordance with theoretical expectation. The helium data on two-proton correlations 
at small relative momentum have been obtained for the first time. The measured size of the 
emission region r^MS depends on A and the pair momentum. Our estimate of the flucton 
size provides a value of rj = 1.3 ± 0.3 fm. 
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